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We construct a class of quantum critical points with non-mean-field critical exponents via holog- 
raphy. Our approach is phenomenological. Beginning with the D3/D5 system at nonzero density 
and magnetic field which has a chiral phase transition, we simulate the addition of a third control 
parameter. We then identify a line of quantum critical points in the phase diagram of this theory, 
provided that the simulated control parameter has dimension less than two. This line smoothly 
interpolates between a second-order transition with mean-field exponents at zero magnetic field 
to a holographic Berezinskii-Kosterlitz-Thouless transition at larger magnetic fields. The critical 
exponents of these transitions only depend upon the parameters of an emergent infrared theory. 
Moreover, the non-mean-field scaling is destroyed at any nonzero temperature. We discuss how 
generic these transitions are. 



I. INTRODUCTION AND SUMMARY 

The AdS/CFT Correspondence [IH3] describes a class 
of strongly coupled gauge theories in terms of weakly 
coupled gravitational systems. It has proved an ex- 
tremely versatile tool for the study of strong coupling 
phenomena over the last ten years. For example, the 
correspondence has been used to study the physics of 
deconfined plasmas, including transport [4] and energy 
loss [5j [6] . Recently there has been much interest in the 
use of AdS/CFT to realize condensed matter phenomena. 
Much of this work has been dedicated to the study of non- 
Fermi liquids [THS] and holographic superfluids [TOlfTT] in 
the hope of better understanding the phase diagram of 
high-temperature superconductors. Another route to the 
same goal involves the study of quantum critical points 
in strongly interacting theories. These zero temperature 
transitions are interesting in their own right, as they tend 
to govern the physics of large swaths of the phase di- 
agram at nonzero temperature. Indeed, the "strange 
metal" phases observed in high-temperature supercon- 
ductors may originate from a quantum critical point |12j . 

The study of critical phenomena is of central impor- 
tance in the condensed matter community. At any con- 
tinuous phase transition there is an emergent infrared 
fixed point |13j . Of particular interest are transitions 
where the infrared theory is itself an interacting quan- 
tum field theory. These transitions are characterized 
by non-mean-field critical exponents. It would be ex- 
tremely interesting if holography can be used to study 
these transitions or perhaps even transitions beyond the 
Landau-Ginzburg- Wilson paradigm altogether [14] . Un- 
fortunately, most continuous transitions in holographic 
models are second-order with mean-field exponents |15j . 
In fact, the mean-field exponents should be expected 
rather than surprising. They appear because of the large 
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N parameter in these theories, which allows us to study 
them via their holographic duals. In this limit quan- 
tum fluctuations are suppressed in both the field |16j 
and gravitational theories. As a result, examples of non- 
mean-field exponents in large N theories are doubly inter- 
esting. Moreover some justification for their non-mean- 
field behavior should be given in the sense of [17] . 

These questions are not only useful for the condensed 
matter community. The study of phase transitions in 
large N theories necessarily sheds light on the physics of 
non-Abelian gauge theories. A general classification of 
transitions in the phase diagram of such theories is im- 
portant. Such a dictionary will help in our understanding 
of QCD-like gauge theories in (3 + 1) dimensions as well 
as in condensed matter systems in (2 + 1) dimensions. 
Much is already known. For example, the finite temper- 
ature deconfinement transitions of these gauge theories 
with holographic duals are first order |18| and map onto 
Hawking-Page phase transitions in the gravity de- 
scription. Additional transitions have been identified in 
systems with quarks 20 27 - there are meson-melting 
transitions in a thermal or high-density bath [3"T1 - 

3D]. Chiral symmetry is also broken in the system with 
a magnetic field 1 1 1!) and there is a chiral restoration 
transition at large densities that occurs in addition to 
the meson-melting transition. All of these transitions are 
typically first order at finite temperature and low quark 
density but continuous at large density (see [50H52"] for 
the full phase diagrams of the D3/D5 and D3/D7 sys- 
tems with magnetic field which we will study here. For 
thermodynamics see [531154] ). 

Flavored gauge theories are also natural to study from 
the critical phenomena perspective. Flavor sectors carry 
new symmetries, leading to a richer phase diagram. At 
large N, we also get the expansion parameter Nf/N 
and so the quenched limit is simultaneously rich and 
tractable. The holographic description of these theories 
involves probe D-branes minimizing their worldvolume in 
the dual geometry [2H] . 

The first example of a holographic quantum critical 
point separating two nonzero density phases was ob- 
tained in the D3/D7 system [501 [5T]. The dual field the- 



ory is simply strongly coupled Af = 4 super Yang-Mills 
(SYM) coupled to a small number of massless fundamen- 
tal hypcrmultiplets. The chiral transition in this theory 
is triggered by large magnetic fields and is second-order 
with mean-field critical exponents. In fact, there are re- 
lated chiral transitions at nonzero density and magnetic 
field for almost all of the supersymmetric probe brane 
systems [55J . Most of these are second-order transitions 
with mean-field exponents, but a handful are not. 

In particular, the first example of a non-second-order, 
non-mean-field transition in holography was identified in 
the D3/D5 system [56 . The dual field theory is the 
same as in the D3/D7 setup, but the flavor fields are 
confined to a (2 + 1) dimensional defect. The chiral tran- 
sition in this theory exhibits exponential scaling and so 
is reminiscent of the celebrated Berezinskii-Kosterlitz- 
Thouless (BKT) transition [57H59] . This is the first 
known instance of exponential scaling at zero temper- 
ature in (2 + 1) dimensions. Indeed, one of us helped 
term this new transition a holographic BKT transition, 
as it occurs in a different context than the original BKT 
transition. Since then holographic BKT transitions have 
also been found in the context of extremal asymptoti- 
cally AdS4 dyonic black holes [50] as well as in two other 
probe brane setups, namely flavored ABJM theory and 
flavored (1,1) little string theory [551151]. They have also 
been identified in noncritical string setups in [62 . 

The existence and properties of the BKT transition are 
intimately related to the Colcman-Mermin- Wagner the- 
orem [53Tf55j . Recall that transitions of the BKT type 
are between disordered and guasi-ordered phases in two 
dimensions. In the quasi-ordered phase, two-point func- 
tions of symmetry-breaking operators have polynomial 
falloff at long distances while the correlation length in 
the disordered phases scales as exp(c/ yJT — T c ) near the 
critical temperature T c [58] . Holographic BKT transi- 
tions are novel in that they exhibit exponential scaling 
in an ordered phase. Correspondingly, their existence 
is entirely unrelated to the long-distance restoration of 
continuous symmetry in two dimensions. On the gravity 
side, holographic BKT transitions occur due to the vi- 
olation of the Breitcnlohner-Freedman (BF) 66 bound 
in the infrared region of the dual geometry by a scalar 
field dual to an order parameter. In the field theory, this 
amounts to considering a theory with an emergent CFT 
in the infrared. The transition is driven by taking an 
operator dimension in the emergent theory into the com- 
plex plane. On general grounds presented in |67j this 
violation was expected to produce BKT scaling, but the 
D3/D5 system was the first example of a setting where 
the BF bound is violated controllably. 

In all known examples of holographic BKT transi- 
tions the dual geometry has an effective infrared AdS2 
region. In the case of extremal asymptotically AdS4 dy- 
onic black holes, the near-horizon geometry is of the form 
AdS2 xK 2 |. For the probe brane systems, there is no 
physical AdS2 region at the bottom of the geometry, but 
worldvolume fields obey the equations of motion for fields 



in an AdS2-like region there [55] • For the dyonic black 
holes the emergent CFT is important: it governs much 
of the low-frequency form of correlation functions at zero 
and small temperatures [9]. In this work, the emergent 
CFT will also be of critical importance. 

With all of this background in mind, there is an impor- 
tant open question: can we use the holographic technique 
to study second-order quantum critical points with non- 
mean-field exponents? If so, can we explain the expo- 
nents at large N? In this work we answer both questions 
in the affirmative 2 . 

Our philosophy is to obtain the most general quan- 
tum critical point in a probe brane system (in the strict 
N — » oo limit). We do this by considering a theory with 
three relevant control parameters at zero temperature, 
one that tends to preserve a symmetry (nonzero den- 
sity) and the other two to break it (one of them is a 
magnetic field). Ideally to explore this one should use 
explicit examples of such a theory. However, the gravi- 
tational description of the probe branes is rather restric- 
tive since only a very small number of operators have 
their dimensions protected and appear as modes in the 
Dirac-Born-Infcld (DBI) action of the brane - the rest 
are stringy modes and have very large dimension. We 
therefore take a phenomenological approach in this pa- 
per and begin with the D3/D5 system at nonzero density 
and magnetic field. We simply include our third control 
parameter, O, into the brane action by hand in a natural 
fashion. It becomes a magnetic field when its dimension 
is two and otherwise we tune its dimension. As long as 
O has dimension less than two 3 , the resulting phase dia- 
gram of the theory is quite rich; we plot it in Fig. [T] At 
fixed density and magnetic field, we find a chiral quan- 
tum critical point as we vary O. Varying the magnetic 
field leads to a line of second-order transitions that con- 
nects to a line of holographic BKT transitions at critical 
magnetic fields. 

In general the exact phase diagram differs depending 
on our exact choice for the dimension of O. However, 
the critical exponents of the chiral transition do not. We 
refer the reader to Figs, [I] and [2] which illustrate these 
points. This "universality" is one of our central results. 
In fact, the critical exponents only depend upon the di- 
mension of the scalar field dual to the order parameter 
in the effective AdS2 region, Air. For example, near the 
transition the order parameter scales as 4> ~ (O — O c )°, 



1 At nonzero temperature, there are two known classes for these 
transitions. In the first, there is a curvature singularity at the 
transition 69 , 70] . The second class is phenomenological, where 
the non-mean-field exponents arise from an effective bulk ac- 
tion 1711 whose terms presumably originate from 1/JV corrections. 

2 As we were finishing this work, the authors of |72| released a 
paper that studies holographic superfluids in the presence of 
double-trace deformations. They obtain quantum critical points 
with non mean-field scalings that precisely match our results. 

3 When O has a larger dimension than the density, chiral symmetry 
is broken for any nonzero O as in 1411 
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FIG. 1: The zero temperature phase diagram of our phe- 
nomenologically deformed D3/D5 system at nonzero (fixed) 
density, magnetic field, and simulated control parameter O of 
dimension A. The shaded region indicates the chirally sym- 
metric phase and the white region the broken phase. For 
A < 2, there is a line of holographic BKT transitions at 
the critical magnetic field B c — There is also a line 

of second-order quantum critical points triggered by O that 
connects to the line of BKT transitions. The position of the 
line depends on the precise manner we introduce O. Indeed, 
the left line corresponds to the choice A = 5/4 and the right 
to A = 1. However, the critical exponents are only a function 
of the magnetic field and density (Figj2|. The static critical 
exponent (3 takes the mean-field value 1/2 for the solid por- 
tion of the lines, but takes on a non mean-field value for the 
dotted section as in Eq. The line of critical points then 
interpolates between a second-order mean-field transition and 
a holographic BKT transition. 
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In the D3/D5 system, Air depends on the ratio of mag- 
netic field to density ( 14 ) , 



1 



A m = 



d 2 -7B 2 
d 2 + B 2 



where B is the magnetic field and d is proportional to 
the density. This dimension goes to unity at zero mag- 
netic field and to 1/2 at the holographic BKT transi- 
tion, B c = d/y/7 (from where it then enters the com- 
plex plane). This is a nice result: our line of transitions 
not only exhibits non-mean-field exponents, but it also 
continuously connects a transition with mean-field expo- 
nents at Air = 1 to a holographic BKT transition at 



Air = 1/2. See also ((62j). 

To understand this interpolation we describe the effec- 
tive potential of the theory near the line of second-order 
transitions. We do this numerically, finding that both 



FIG. 2: The critical exponent /3 in the deformed D3/D5 sys- 
tem as a function of magnetic field at A = 1 and A = 5/4. For 
small magnetic fields - such that Air > 3/4 - the exponent j3 
assumes a mean-field value, while for larger magnetic fields it 
does not. The exponents for A = 1 (blue dots) and A = 5/4 
(red dots) match each other and our prediction (dotted line), 
Eq. (621, within our numerical accuracy. 



the order parameter and free energy in the broken phase 
follow from a modified Landau-Ginzburg model with a 
potential of the form 



V e f}{4>) = a 2 {O c - <J)4> 2 + cv 4 </> 4 + am(f> 



l/(l-Ai R ) 



(2) 



where the cv's are positive, dimensionless couplings and 
Air is bigger than 1/2 and less than 1. When O is tuned 
past its critical value O c , the (f> 4 dominates for Air > 3/4 
and the last term dominates for Air < 3/4. In the sec- 
ond regime, the static critical exponent (3 takes on a non 
mean-field value. The non-analytic term in the poten- 
tial has a natural form given that there is an emergent 
(0 + l)-dimcnsional infrared theory under which the con- 
densate has dimension 1 — Air. The term is just that 
required on dimensional grounds in the infrared theory. 



We justify this further in Sec. IV 



Our thermodynamic results are obtained numerically. 
However, we do obtain some analytic results for fluctua- 
tions in the symmetric phase of the theory. Following [9] , 
we perform a matching computation to obtain the low- 
energy limit of the two-point function of the order pa- 
rameter. The result for the brane system is essentially 
the same there. Assuming some basic analyticity con- 
straints, we compute both the correlation length and the 
dynamical critical exponent near our transitions. The 
former diverges with mean-field scaling while the latter 
assumes a non mean-field value for all nonzero values of 
the magnetic field. As with /?, it only depends on Air . 

What happens at nonzero temperature? In the pre- 
vious examples with both probe branes and holographic 
superfluids, the transitions are universally second-order 
with mean-field exponents [56] . Even BKT scaling is de- 
stroyed at any nonzero temperature. For this reason, we 
also expect the non-mean-field scaling of our transitions 
to be lost away from zero temperature. Indeed we find 



that this is the case below. 

Our general conclusion then is that a conformal the- 
ory, perturbed by three control parameters O s ,0\ 2 with 
dimensions A s = Ai > A2 may lead to a phase dia- 
gram qualitatively similar to that represented in Fig. [T] 
In this general picture we require that O s tends to re- 
store a symmetry and 12 to break it. The constraint 
that A s = Ai amounts to the freedom to change the 
emergent infrared theory in these systems. Meanwhile, 
the Ai > A2 condition allows us to trigger an ordinary 
second-order transition without altering the infrared the- 
ory Since our results crucially depend upon an emergent 
CFT, we expect that they extend beyond probe branes 
to all systems with such an emergent theory. 

We further test this picture phenomenologically by 
performing a similar analysis of the D3/D7 system. In 
that theory, the magnetic field has dimension two and 
the baryon density dimension three. We show that if a 
chiral symmetry-breaking dimension three operator, O, 
is also introduced, then the theory realizes a holographic 
BKT transition at large O. As in the D3/D5 system, in- 
troducing a control parameter with the same dimension 
as the density leads to an emergent theory where the di- 
mension of a scalar operator depends on the density and 
O. Moreover, at intermediate O the magnetic field can be 
used to trigger a second-order non-mean-field transition. 
While we have not computed it numerically, we expect 
that the phase diagram qualitatively matches our results 
for the D3/D5 system. 

The outline of this work follows. In Sec. |TT] we re- 
view the D3/D5 and D3/D7 systems at nonzero density 
and magnetic field [50—52, 56 . We go on to present our 
phenomenological models and our holographic regular- 
ization. The bulk of our results are presented in Sec. |III[ 
beginning with the thermodynamics of the D3/D5 sys- 
tem at zero and nonzero temperature. Next, we study 
fluctuations using a matching procedure. In Sec |IV| we 
present our effective theory of the transition and criti- 
cally test it. We apply our analysis briefly to the D3/D7 
system without numerics in SecJVj Finally, we discuss 
our results in Sec. IVIl 



II. HOLOGRAPHIC SETUPS 

A. The D3/D5 system 

Strongly coupled SU(N) W=4 super Yang-Mills 
(SYM) theory at large N and zero temperature is dual to 
type IIB supergravity on AdSs x 5 5 (with N D3 branes 
at its core). The geometry can be written as 

2 jd2 

ds2 = R?^' 1 + w 2 ^ + + ^ + L2 ^2): (3) 

where w 2 = p 2 +L 2 1 dil 2 , dfl 2 , are the metrics for two unit 
two-spheres, and R A = 4irg s Na 2 . In these coordinates 
the Poincare horizon of AdS is located at p = L = and 



the boundary as p 2 + L 2 — >• 00. 

We now add Nf flavor hypermultiplets to the gauge 
theory along a (2 + 1) dimensional defect by placing a 
probe D5 brane in this geometry. The probe limit corre- 
sponds to the quenched limit of the gauge theory. The 
D5 probes are described by their DBI action 

Sdbi = -N f T 5 J d 6 ^-dct(P[G] ab + F a b) , (4) 

where a, b = 0,..,5 are world volume indices, P[G] a b 
is the pullback of the metric to the brane, and F is 
the field strength for the diagonal U(l) gauge field liv- 
ing on the D5 worldvolume. The field theory has a 
50(3)1 x 50(3)2 x U(1)b global symmetry, where the 
two 50(3)'s are chiral R-symmctrics and the U(1)b is 
a baryon number symmetry which only rotates the fla- 
vor fields. The baryon symmetry current is dual to the 
U(l) gauge field on the brane while the chiral symmetries 
correspond to the rotational symmetry of the two two- 
spheres. We now consider an ansatz wherein our brane 
embeddings are translationally invariant in the wrapped 
field theory directions x° — x 2 while wrapping the first 
two-sphere and the "radial coordinate" p. We also con- 
sider the theory with no source for the baryon current, 
that is with F = 0. There are then a family of embed- 
dings L — m and x 3 ,92,4>2 constant (where #2>^>2 are 
coordinates on the unwrapped two-sphere). These are 
supersymmetric embeddings that correspond to the the- 
ory with a hypermultiplet of mass m. The second 50(3) 
chiral symmetry is then explicitly broken by a hyper mass 
and spontaneously broken by a vev for the correspond- 
ing operator tptp (plus operators related by supersymme- 
try) |23) . It is this chiral symmetry that is spontaneously 
broken in our chiral transitions. 

We now extend our ansatz to include nonzero baryon 
density and magnetic field. To do this, we need to let the 
embedding function L depend on p as L — L(p) as well 
as turn on a non-trivial field strength [Kl EI] , 

F = A' (p)dp A dx° + B dx 1 A dx 2 , (5) 

where the field Aq will determine both the chemical po- 
tential and density of baryon charge [371 EH EH [74] and 
B is the magnetic field. The radial electric field A' is 
sourced by charge at the bottom of AdS, and so the 
brane ends there. This amounts to the boundary con- 
dition L(0) = 0. For this ansatz we can consistently 
neglect the Wess-Zumino pieces of the brane action and 
write 

5,5 = - N f T 5 R 6 vo\[S 2 } volpR 2 ' 1 ] 

" v ' 

x J d pp 2 ^l + L' 2 -A' 2 ^l + ^ (6) 

where we have defined w 2 = p 2 + L 2 as well as rescaled 
x and p by powers of R. The normalization is given 



by TV = y/XN/2ir 3 where A = Ang^ M N is the 't Hooft 
coupling of the SYM theory. From here onward, we will 
refer to the action density S^/volpR 2,1 ], describing it with 
the same notation S5. 

Notably, the action only depends on Ag through its 
radial derivative. Thus there is a conserved quantity, 



d = 



SS* 



1 + 



B 2 



SA' (p,x) ^1 + L' 2 - Aj 



(7) 



In fact, d is the baryon density in the dual theory. Solving 
for A' in terms of a rescaled density d = d/J\f, we find 



A' 2 - - 



d 2 {l + U 



d 2 + p 4 (1 



B 2 > 



(8) 



We obtain the brane action at fixed density by substitut- 
ing this result into the action Eq. [6] and Legendre trans- 
forming with respect to A' . The result is 



S 5 = -Af J d P Vl + L' 2 Jd 2 + P i (l 



B 2 

w 4 



(9) 



Field configurations L{p) that extremize this action cor- 
respond to field theory ensembles that extremize the ef- 
fective potential of the theory in the canonical ensemble. 
In general, these configurations can only be obtained nu- 
merically. However there is an exact solution to the equa- 
tion of motion for all d and B, L = 0, which corresponds 
to the dual theory with zero hyper mass and zero conden- 
sate. This solution corresponds to the chirally symmetric 
phase of the theory. 

We continue by reviewing the origin of the chiral BKT 
transition in this system. The onset of the transition 
can be understood by studying the stability of the sym- 
metric embedding. Small fluctuations around L = are 
described by the quadratic piece of Eq. ^ , 



C 5 



d 2 + B 2 + p 4 L' 2 



MB 2 L 2 



(10) 



d 2 + B 2 



This Lagrangian has two distinct limits 4 . At large p ^> 
\/B, yd, the field L/p fluctuates as a stable m 2 = —2 
scalar field in AdS4. However, at small p <C VB, vQ, 



4 In our analysis we use the results for a scalar in AdS p +i: The 
solution of the equation of motion is 



L 

P 



(11) 



(12) 



and the Breitenlohner-Freedman (BF) bound |66] is given by 
-p 2 /4 



L/p fluctuates as a m 2 = —2B 2 /(d 2 + B 2 ) scalar in AdS2- 
Thus for d/B > V7 the field is stable but for d/B < sfi 
the mass drops below the Breitenlohner-Freedman (BF) 
bound 66 in AdS2, ™|f = — 1/4- There is therefore a 



chiral transition at the critical filling fraction 



B, 



7XN f N 



(13) 



Further analysis reveals that the order parameter scales 
exponentially at smaller densities, so that the transition 
is of the holographic BKT type. This behaviour is the 
result of the violation of the BF bound in the infrared 
region which implies that an infinite number of tachyons 
form at the transition, an extremely unnatural situation 
within the Landau-Ginzburg- Wilson paradigm. Because 
it will be important in the rest of this paper, we also note 
that in the effective AdS2 region L/p is dual to a scalar 
operator in the emergent CFT of dimension 



1 



A 



IR 



d 2 -7B 2 
d 2 +B 2 



(14) 



As usual in AdS / CFT there are two solutions to the equa- 
tion of motion for the scalar in AdS (see Eq. ( 12 )) which 
describe an operator in the field theory and its source. 
The second solution here corresponds to an object of di- 
mension of 1 — Air. We will see below when we discuss 
the effective theory for our transitions that the 3d con- 
densate corresponds to a dimension 1 — Air operator in 
the infrared theory. 

At zero magnetic field we have Air = 1, which de- 
creases to Air = 1/2 at the transition. In the broken 
phase, Air is driven into the complex plane. In this way 
the scaling symmetry of the infrared theory is broken to 
a discrete subgroup (which is broken further to a self- 
similar subset by higher energy physics), which relates 
the various tachyons of the symmetric phase. 



B. The D3/D7 system 

In the same way we can consider strongly coupled 
N = 4 SYM at large N coupled to (3 + 1) dimensional 
fundamental hypermultiplets. In the quenched limit the 
flavor fields are well described by probe D7 branes in the 
AdSs x S 5 geometry. The global symmetry of this theory 
is 50(4) x U(l) x x U(1) B , where the U(l) x is a chiral 
symmetry and the U(1)b is the usual baryon number 
symmetry. This chiral symmetry is explicitly broken by 
a hyper mass and spontaneously broken by a condensate 
of the hyper mass operator. 

On the gravity side, the D7 branes wrap a three-sphere 
rather than a two-sphere. The U(l) x chiral symmetry is 
dual to the 5*0(2) isometry of an K 2 transverse to both 
stacks of D3 and D7 branes. The baryon symmetry cur- 
rent is dual to the diagonal U(l) gauge field on the D7 
branes as before. We then consider a translationally- 



invariant, S'0(4)-preserving ansatz as before with a den- 
sity and magnetic field. For such an ansatz the brane 
action at fixed density is 



S 7 = -A/V / dp y/l + L^A dP 



B 2 

w 1 



(15) 



As above, we can study the onset of the chiral transi- 
tion by studying the stability of small fluctuations around 
the chirally symmetric embedding L = 0. These are de- 
scribed by the quadratic part of Eq. ( 15 ), 



96L' 2 



N 7 B 2 L 2 



d 2 + p 2 B 2 



(16) 

As with the D3/D5 system, this Lagrangian has two dis- 
tinct limits. For p ^> d 1 ^ 3 , yB, L/p fluctuates as a stable 
m 2 = — 3 scalar field in AdSs. On the other hand, for 
p <C d 1 / 3 , VB L/p fluctuates as a stable massless scalar 
in AdS2- As originally pointed out in [55] . there is an 
emergent CFT in this theory as well. Moreover, L/p is 
dual to a scalar operator in the infrared theory of dimen- 
sion Air = 1. 

There is no holographic BKT transition in this system. 
Rather, the chiral transition is second-order with mean- 
field exponents. A single tachyon forms at the transi- 
tion, which is effectively modelled by a Landau-Ginzburg 
model with a quartic potential. Later, we will see that 
the mean-field exponents are crucially related to the fact 
that Air = 1. For the majority of this paper we will 
work in the D3/D5 system that does have a BKT tran- 
sition but we will return at the end to produce similar 
phenomena in a phenomenological deformed version of 
this D3/D7 system. 



C. Phenomenological models 

We seek to extend the brane systems above by turning 
on a third control parameter. For computational sim- 
plicity, we seek to deform our setups in such a way that 
the brane action depends only upon a single worldvolume 
field and a number of constants of the motion. For the 
D3/D5 system, there are several candidates. The first 
is an electric field along the brane [75] and the second a 
flux on the wrapped two-sphere [76] . Neither deforma- 
tion breaks chiral symmetry at zero magnetic field and 
zero hyper mass, but the electric field may yet lead to in- 
teresting results. Other deformations involve additional 
worldvolume fields whose equations of motion are not in- 
tegrable. 

In favor of solving a more complicated brane problem 
with at least two worldvolume fields, we elect to take 
a phenomenological approach. We will simulate a chi- 
ral symmetry-breaking control parameter whose dimen- 
sion we dial. At first, this approach may seem cavalier: 
in contrast with the "bottom-up" holographic superfluid 



and non-Fermi liquid analyses, there are many different 
ways that control parameters emerge in a probe brane 
action. There are few a priori reasons to believe that 
phenomenology will accurately predict features of tran- 
sitions in consistent "top-down" brane setups with three 
control parameters. 

The best justification for our method comes ex post 
facto. Ultimately, we find that the critical exponents we 
measure do not depend upon the details of our simulated 
deformation. This result is crucial and we will return to 
it extensively later. For now, we will simply describe our 
phenomenological choice. We simulate a control parame- 
ter O of dimension A (taken to be relevant) in the D3/D5 
system by considering a modified brane action 



S 5 = -JV J dp^l + L' 2 ^ 2 + p 4 (l + ^ + ° 2 

(17) 

Note that when O has dimension two, it is effectively a 
magnetic field. 

Earlier we studied both the onset of the chiral tran- 
sition as well as the emergent infrared theory by study- 
ing small fluctuations around the symmetric L = em- 
bedding. These are now described by the quadratic La- 
grangian, 
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(18) 



This system has two different infrared limits depending 
upon the value of A. For A > 2, the new control pa- 
rameter dominates the infrared, so that L/p fluctuates 
as an unstable scalar there. Then for any nonzero O the 
symmetric embedding is unstable and the stable phase 
will break chiral symmetry. For A < 2 the magnetic 
field and density together dominate the infrared. As be- 
fore, at small p <C VB, Vd, the field L/p fluctuates as a 
m 2 = —2B 2 /(d 2 + B 2 ) scalar in an effective AdS2 region 
at the bottom of the brane. Finally, the introduction of 
O tends to break chiral symmetry. We see this by study- 
ing small fluctuations at nonzero O but vanishing density 
and magnetic field. In this limit, the radial equation of 
motion for the field L/p at small p is that of a m 2 = A — 3 
scalar in AdS4_A- Provided that O is not marginal with 
A = 3, L/p fluctuates unstably in the infrared and so the 
true ground state breaks chiral symmetry as claimed. 

We can make a similar phenomenological deformation 
to introduce an operator O of arbitrary dimension into 
the D3/D7 system. Just as in the D3/D5 system, we 
introduce a simulated control parameter into the brane 



action as 



define a regularized action by 



SV = -A/7 J dp y/l + L' 2 Jd 2 + p 6 ^1 



B 2 O 2 



(19) 

We will discuss the physics of this model in the later 
Sec|Vl 



D. Holographic regular ization 



The phenomenological brane actions Eq. ( 1 7 ) , ( 19 ) con- 
tain a number of near-boundary divergences [77\. In a 
genuine "top-down" construction, these correspond to 
ultraviolet divergences of the dual theory. In the bulk, 
they can be diffeomorphism-invariantly regulated by in- 
troducting a near-boundary cutoff slice and adding local 
counterterms on the slice. This process is known as holo- 
graphic renormalization and is crucial: once we have an 
appropriately renormalized bulk action, we may sensibly 
take derivatives to obtain correlation functions. 

In our phenomenological constructions, however, we 
cannot be sure that we diffeomorphism-invariantly regu- 
late the bulk theory. Rather we choose to regularize our 
theories in a manner inspired by holographic renormaliza- 
tion. We introduce a cutoff slice, add counterterms, and 
then take the cutoff to infinity. To illustrate the idea, we 
consider the modified D3/D5 system of Eq. (17 1, where 
O has dimension A — 1. A general solution L[p) has the 
near-boundary solution 



L(p) = m + J2-^-, 

71=1 H 



(20) 



where the L n for n > 1 are recursively determined by 
m and L\. The parameter m is simply the hyper mass. 
Then the brane action, integrated up to a cutoff p = 
A, evaluated on such a solution has the near-boundary 
divergence structure 
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N / d 3 x 
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2 A 



finite. 



(21) 



The exact divergences depend upon A; for A > 3/4 
there is only a single counterterm required, while for 
A = 3/(2n) for n a positive integer there is also a loga- 
rithmic divergence. These logarithmic divergences corre- 
spond to Weyl anomalies of the dual theory [78] . In this 
case of A = 1, we add simple counterterms on the cutoff 
slice, 
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(22) 



where 7 is the induced metric on the slice. When there 
are logarithmic divergences, we subtract them with coun- 
terterms of the form ^—^0° / (p 3 log p) . In general, we 



^.reg = lim 

A— too 



£>5,A + Si 



CT 



(23) 



We now define correlation functions of the dual theory 
through functional derivatives. 

For A > 1/2 there are no divergences that depend on 
771 or L\ . Consequently in this region the one-point func- 
tion of the operator dual to L, the hyper mass operator 
Ol, is simply 



<t> - (Ol) 



SS: 
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(24) 



On the other hand, for A < 1/2 there are additional 
contributions to (O y ) that arise from extra L-dependent 
counterterms. These counterterms are proportional to 
L(p) 2 and so lead to a contribution to the condensate 
proportional to 777. These contributions vanish in the 
chiral limit we consider. 

A similar analysis can be performed for the D3/D7 
system. Since we do not use the results in this work, we 
simply quote the highlights. For the choice A = 3 that 
we consider in this work, there are only two divergences: 
the first corresponds to the infinite volume on the cutoff 
slice and the second to a logarithmic divergence propor- 
tional to Fu, v F^ v evaulated on the slice. The regulariza- 
tion is then identical to the holographic renormalization 
performed in |51j . Notably, the logarithmic divergence 
corresponds to a Weyl anomaly of the dual theory: the 
trace of the stress tensor of the dual theory is propor- 
tional to F 2 . 



III. RESULTS 

In this section we report the bulk of our numerical re- 
sults. We begin by studying the zero temperature transi- 
tions of our deformed D3/D5 system. As a first numeri- 
cal check we reproduce the BKT transition in the D3/D5 
system with a magnetic field and density. Next, we study 
this system with the control parameter O present but at 
zero magnetic field and show that it triggers a mean-field 
chiral transition. Next, we study the zero temperature 
phase diagram in the (O, B) plane and identify the line of 
quantum critical points. We go on to study the nonzero 
temperature transitions, for which the non mean-field 
scaling is destroyed. Finally we compute the low-energy 
behavior of the two-point function of the order parameter 
in the symmetric phase. 



A. Zero temperature transitions 

1. BKT transition 

We begin by studying the D3/D5 system with mag- 
netic field and density to reproduce the BKT transi- 



tion [SS] . To do so we must solve the equation of motion 
for L(p) that follows from varying Eq. (17) with O zero. 



This equation can be solved numerically. We do so with 
the shooting method, generally shooting from large p. 

Recall that the radial electric field on the brane is 
sourced by charge at the bottom of the geometry. This 
is equivalent to setting the infrared boundary condition 
L(0) = 0. At small p there is then a series solution for L 



HP) = 7o + liP + 7n 



(25) 



where the higher 7„ are recursively determined by 70 and 
7i- 

Near the AdS4 boundary we impose the boundary con- 
dition that our dual flavor is massless. This amounts to 
choosing the leading term in the near-boundary solution 
Eq. ( 20 1 to vanish. For our shooting we use the large p 



series solution as initial data (having computed the first 
dozen or so L n s) upon which we numerically integrate 
the corresponding solution to small p. We then match 
this solution onto the small p solution. By dialing the 
field theory condensate we shoot for solutions that ex- 
tend to the bottom of AdS with 70 = 0. It is numerically 
difficult to match onto a small p series solution. A sim- 
ilar problem emerges in the small p embeddings in the 
flavored little string theory studied in [55]: in each case, 
there is a non-analyticity in the equation of motion for 
L at p = 0. For the D3/D5 system, embeddings conse- 
quently "spike" to p = 0, infinitely so at infinitesimally 
small p. Nonetheless, with care, we have managed to suc- 
cessfully shoot to the infrared boundary condition with 
high accuracy. 

Generically with large B the solutions bend off the L = 
axis. The chemical potential then forces the solution 
to spike to the origin at p = L = 0. At some critical 
value of B the embeddings smoothly transition to the 
L = embedding. These embeddings are shown in Fig 
[3] for varying B/d. We also plot the value of the quark 
condensate, c, across the embedding to show the BKT 
exponential scaling. 

Note that at yet larger B there is a second transition 
to a phase with d — and stable mesons described by 
embeddings that curve off the axis but do not spike to 
the origin. That transition is described for this theory in 
[52) but we will not explore it further in this paper. 



2. The phenomenological operator O 

Next we study the deformed D3/D5 system at zero 
magnetic field. As we discussed in Sec. |II C when the 
dimension of O is greater than two we find that the chi- 
rally broken phase is preferred for all values of the den- 
sity. This matches the prediction above that the L = 
embedding is unstable for all d. 

When the dimension of O is less than two, we expect 
large O at fixed density to trigger a chiral transition. We 




(a)The embedding L of a D5 brane in the D3 
geometry for various B/d showing the BKT 
transition. 
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(b)A plot of the quark condensate c versus B 
across the D3/D5 BKT transition. 



FIG. 3: The BKT transition in the D3/D5 system with quark 
density and magnetic field present. 



seek to locate this transition for many different dimen- 
sions A and to measure the associated critical exponents. 
To do so we must solve the equation of motion for L(p) 
that follows from varying Eq. (17) with B = 0. Again 
we use the numerical techniques discussed in Sec. Ill A 1 
above. 



We find that for large O ^> d A / 2 , the solutions cor- 
responding to zero mass bend away from the symmetric 
embedding L = 0. Near the AdS4 boundary they nec- 
essarily asymptote to the symmetric embedding, but at 
small p they "spike" to the bottom of AdS. Thus chiral 
symmetry is indeed broken at large O as expected. 

For fixed A < 2, there is a critical O c where the embed- 
dings smoothly transition to the symmetric embedding. 
This is the location of the chiral transition. We locate 
the transition for many different A, for which we must 
also scan through O. This is somewhat laborious, as we 
have to shoot for each value of O and A. The net result 
is shown in Fig. [4] For A < 2 we identify a second-order 
transition with mean-field exponents as promised. At 
A = 2 there is a holographic BKT transition, as O then 
acts like a magnetic field. We plot the condensate near 
the transition for the particular case A = 1 in Fig. [5] to 
display that mean field behavior. 
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FIG. 4: The zero temperature "phase diagram" of the de- 
formed D3/D5 theory at zero magnetic field as a function of 
the dimension and value of the deformation O. The shaded 
region is the chirally symmetric phase and the white is the 
broken phase. The line of transitions is second-order with 
mean-field exponents, excepting a holographic BKT transi- 
tion at A = 2. 



transitions. The resulting phase diagram as a function of 
O and magnetic field was already plotted in Fig. [T] We 
also plot the critical exponent /3 as a function of magnetic 
field in Fig. [2j Recall that j3 is the scaling of the order 
parameter in the broken phase, c ~ (O — O c )P . 

The combined results are very interesting. First, the 
line of transitions is second-order as expected. Moreover, 
for small magnetic fields such that Air > 3/4, the ex- 
ponent (3 takes on the mean- field value 1/2. Once the 
infrared dimension dips below 3/4, however, the expo- 
nent /3 is no longer 1/2 and moreover is independent of 
the dimension of O. The simplest way to interpret this 
result is that the effective potential of the theory near the 
transition has the usual quartic term as well as a second 
term that depends upon Air but not A. We will show 



how this occurs explicitly in Sec. IV where we construct 
a modified Landau-Ginzburg model for the transition. 



B. Nonzero temperature thermodynamics 




It is interesting to also study the behaviour of our 
model at nonzero temperature. At nonzero temperature, 
the AT = 4 SYM theory is holographically described by 
IIB supergravity on an AdSs black brane geometry (with 
N hot D3 branes at its core). The geometry can be writ- 
ten as 



ds 2 



where 



fH = i- 



%{-f{w){dx»f + dx 2 ) + jj^d^ + tfdn* 
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FIG. 5: The condensate in the deformed D3/D5 theory at 
zero magnetic field and A = 1. The solid line is numerical 
data and the dotted a fit. Near the transition the condensate 
scales with a mean- field exponent, c ~ y/O — O c . 



3. The general case 

Now we turn on the magnetic field and O together. 
This will tune the dimension A/^ of the operator dual to 



L/p in the emergent theory as we discussed in Sec. II C 



As above, we can only solve the equation of motion for 
L numerically. Our procedure is essentially the same as 
the one described above in Sec. IHI A II 

We studied two different dimensions A for O in great 
detail, A = 1 and A = 5/4. The choice of a non- integer 
dimension explcitly shows the independence of our re- 
sults on that dimension. In each case, we studied the 
chiral transition at many different magnetic fields and 
consequently at many values of O as well. For each such 
choice of O and B we must shoot to find the correct vac- 
uum. As expected, we find a line of second-order chiral 



and we define Wh = ttT/R 2 with T the temperature of 
both the field and gravitational theories. This coordinate 
system is related to the one we employ at zero tempera- 
ture by L — wsm9, p = wcosO. Both 9 and L are dual 
to the hyper mass operator. We change coordinates sim- 
ply because we have found the numerics easier for this 
analysis. 

As before, we embed Nt D5 branes in this geome- 
try. We consider embeddings that are translationally- 
invariant in the wrapped x° — x 2 directions, wrap the first 
two-sphere and w, and possess no angular momentum on 
either two-sphere. The embedding is parametrized by 
the worldvolume field 6 — 8(w). After adding a charge 
density and magnetic field, the D5 action at fixed density 



S b = -M J dw^/l + fw 2 6' 2 ] jd 2 + u) 4 cos 4 6 (\ + 
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(28) 



Now we add our third control parameter O. There is yet 
further ambiguity in how we phenomenologically intro- 



duce O. We elect to add it in such a way that it again 
becomes a contribution to the magnetic field when the di- 
mension of O approaches two. Our deformed Lagrangian 
is 



c 5 = -A/yi + fw 2 e ,2 Jd 2 



B 2 O 2 



„2A 



(29) 

At zero temprature the holographic BKT Transition 
was triggered by driving the mass of L in the effective 
AdS2 region below the BF bound. At any nonzero tem- 
perature this exponential scaling is lost [521 [53] ■ The 
infrared region becomes an AdS2-like space with a black 
hole, which has a Rindler near- horizon limit. Driving the 
mass below the BF bound then corresponds to a UV in- 
stability (but not an IR instability) from the point of view 
of the infrared theory. This instability is tamed by the 
ultraviolet completion to AdS4 physics. The absence of 
such an IR instability presumably leads to the resulting 
second-order mean-field transition observed at extremely 
small temperatures in the D3/D5 system. 

By the same logic we expect the non mean-field scal- 
ing of our second-order transitions to be destroyed at 
any nonzero temperature. Indeed, we find this result 
numerically. To do this, we extremize the modified ac- 
tion Eq. ( 29 1 and regularize the bulk action in such a 



way that we measure the field theory condensate from 
our embeddings. As at zero temperature, we employ a 
shooting technique. This time we elect to shoot from the 
infrared near the black brane horizon. The charge on the 
brane indicates that the embedding extends down to the 
horizon. Our infrared boundary condition is then simply 
that the embedding is regular there. 

We plot the condensate near the nonzero temperature 
transition with small temperature irT = l0~ 5 \/d, mag- 
netic field B = 0.98B C , and the choice A = 1 in Fig.[6] At 
zero temperature and this magnetic field the condensate 
scales with an exponent f3 = 2.18, noticeably different 
from the mean-field value obtained at the small temper- 
ature shown in the figure. 



C. Fluctuations 

We now move on to consider dynamics. In particular 
we will obtain the structure of the retarded two-point 
function of the order parameter at low frequency. This 
computation essentially mimics the scaling and match- 
ing methods employed in [S], so we will only quote the 
highlights. 

We begin by considering a fluctuation of the world- 
volume field 8 around the symmetric embedding 8 = 0. 
There is an infinite tower of AdS4 modes correspond- 
ing to the Kaluza-Klein harmonics of 8 reduced on the 
wrapped two-sphere; we only consider the s-wave, as this 
is the lightest mode in the tower. The chiral transition 
destabilizes it. 




FIG. 6: The condensate in the deformed D3/D5 system at 

zero (solid) and small (dashed) temperature ttT = 10 
large magnetic field B = 0.98B C , and the choice A — 1. The 
non mean-field scaling at zero temperature is destroyed even 
at this small temperature. The nonzero temperature conden- 
sate asymptotes to the zero temperature value far away from 
the transition. 



After a straightforward, if frustrating, computation we 
obtain the full Lagrangian for a time and spatially depen- 
dent embedding 8 = 8(x a ,x l ,w). Denoting 



8' = d w 6, 8 = d x o8, (W8) 2 = (di0)(diO), 



(30) 



where i = 1,2, we find 
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(31) 



Now we must choose how to implement our deformation 
for the system with a spatially-varying 8. As before, we 
make the choice that O becomes a magnetic field for A = 
2. This amounts to taking 
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(32) 



The two-point function of the condensate may be com- 
puted by solving the bulk action to second order in the 
variation. That is, by solving the linearized problem 
around 8 = 0, 



V8 = 0, 



(33) 



where T> is & nasty second-order differential operator. We 
also Fourier transform in the x* 1 directions and impose 
the incoming boundary condition at the horizon. The 



resulting solution will have a near-boundary expansion 



regions 
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(34) 



Moreover, the two-point function of the condensate is 
computed by varying the regularized bulk action twice 
with respect to 6±, 



G(u, k) = (0(lo, k)0(~co, -k)) = K 



i(ui,k)' 



(35) 



for K a positive constant. Solving for G at low energies is 
a bit tricky as the correct infrared behavior of 9 depends 
on uj at leading order. We therefore solve 9 in the infrared 
region and match it to the outer region. 



In order to solve for 9 near the bottom of the brane we 
employ a scaling limit. Consider 



W = — , t = A T, 



(36) 



in the A — > limit with £, r finite. At zero temperature, 
the equation of motion of 9 becomes that of a m 2 = 
-2B 2 /(d 2 + B 2 ) scalar field in AdS 2 , 

& e +(u 2 + J I3 ^-±)d = 0. (37) 
? V <P + B 2 e) 

At nonzero temperature, however, we supplement the 
scaling limit Eq. ( 36 ) with 



Wf, 



£o finite. 



The equation of motion in the infrared is then 
P 2B 2 1 



where 
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= 0, 



(38) 



(39) 



(40) 



At zero temperature, the scaling limit Eq. ( 36 1 amounts 
to the uj <C yJ~B limit. At nonzero temperature 
Eq. Q is the uj,T < Vd, VB limit with oj ~ T. 
Notably we can solve for 9 in this region at both zero 
and nonzero temperature. Even more importantly, these 
equations of motion are (i.) independent of both O and 
its dimension, A and (ii.) those of a scalar in either AdS 2 
or an AdS 2 space with a black hole. 



Unfortunately the scaling limits Eq. ( 36 ) , ( 38 1 do not 
give rise to a systematic matching program. As noted 
in [5], a proper matching divides the w axis into two 



Inner: w = —, for £ G (e, oo) 
Outer: — < w, 



in the limits 



0, £ = finite, e -> 0, 



0. 



(41) 
(42) 

(43) 



Small uj perturbations can be treated systematically in 
each region, employing £ as the coordinate in the inner 
one and r for the outer. The result has the form 



Inner: = <?f (£) + 0^(0 + . . . 

Outer: 9 {r) = 9 < g ) (w) + uj9^ ] (w) + 



(44) 
(45) 



The domain of these solutions overlaps in the region de- 
fined by £ — > with w = u/£ —> 0; the full solution is 
obtained by matching 9j and 9o there. 

Now we solve for 9 in the inner region. The leading or- 
der equation of motion for is identical to the one we 
found after the scaling limit, Eq. (39). Near the bound- 
ary of the AdS 2 region (that is, £—►()), the leading order 
term in 9j can be expanded as 



(46) 



where y>±(£) are the non-normalizable/normalizable so- 
lutions to Eq. (391 and Gir{uj) is the retarded Green's 



function of the operator dual to 9 in the infrared the- 
ory. It takes on two vastly different forms depending on 
whether we are at exactly zero or nonzero temperature. 
For the first, it is 



g Am (u>) ex M 2A »-\ 

while at nonzero temperature it is 

C/AmH cx (ia;)T 2AlR_1 . 
The precise form of Q can be found in [9]. 



(47) 



(48) 



The bottom of the outer region corresponds to the 
near-boundary region on the infrared AdS 2 . The solu- 
tion to 9o therefore has the same functional form there, 
and so we can choose a basis where the linearly inde- 
pendent solutions for 9o match precisely to <p± in the 
infrared. That is, 

9^\w) = rf+\w) + GaM V ^(w), (49) 

where r]± ^ (w) is our (zeroth-order) basis in the outer re- 
gion. At higher order in uj the matching can be system- 
atically employed, effectively correcting the basis at each 
order so that 



9 (w) = r)+(w) + g AlR (u)r)-(w) 



(50) 



is satisfied. 

Near the AdS4 boundary the n th order corrections to 
r]± will have an expansion 
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This together with Eq. (35) leads to the desired result, 



namely the form of the retarded two-point function 
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Moreover, by expanding the a's and 6's about fe = 0, we 
find that G assumes the small w, k form 
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as long as Air < 1. For Air > 1, the low-energy limit is 
instead 



G(tu,k) 



g a + giu + g 2 k 2 
7o + fiu + hk 2 



(54) 



On the reasonable assumption that the matching coef- 
ficients, gi and fi, are analytic in O — O c , then the chiral 
transition corresponds to a root in /o. That is, near the 
transition /o is proportional to O c — 0. We then find that 
for Air < 1 the mode that drives the transition obeys a 
zero temperature dispersion relation 



2A m -l 



f 2 k 2 <xO c -Q. 



(55) 



From this relation we simultaneously obtain the dy- 
namical critical exponent z at the transition, 



(56) 
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and the divergence of the correlation length, 
(0(aO0(O)> ~ e-W/e, e~(Oc-0)-", v 



(57) 



Thus the dynamical exponent takes a non mean-field 
value while v takes the mean-field one. At nonzero tem- 
perature, however, the dispersion relation becomes 



fo(iu)T 
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hk 2 = O c -Q, 



(58) 



so that the dynamical critical exponent takes on a mean- 
field value z = 2. As with the condensate, the non mean- 
field scaling is destroyed at any nonzero temperature. 



IV. AN EFFECTIVE MODEL 



there is a rich phase structure. The chiral restoration 
transition is of the holographic BKT type at large B but 
second order at large O with a region in between with 
non mean-field behaviour. Our results provide a rich ar- 
ray of numerical data that we will show can be completely 
matched by a simple effective theory. 

We 5 have been able to guess the form of the effective 
potential. For Air > 3/4, the mean field scaling of static 
exponents is reproduced by the potential 

VeffW = V + a 2 {O c - 0)4> 2 + « 4 4 + O(0 6 ), (59) 

where the couplings on are presumed to be positive and 
Vo is the free energy in the symmetric phase. This gen- 
erates the expectation value <j> ~ \JO — O c .This is just a 
standard Landau-Ginzburg model. 

The crucial ingredient when we move away from mean- 
field scaling is that the gravity dual reveals that the infra- 
red dynamics is governed by a lower dimensional AdS p +i 
theory (for the probe brane systems, p = 1). We have 
also learnt that the order parameter in this low energy 
regime acts as an operator of either dimension A/# or 
p — Ajfj (see Eq. (14) for the D3/D5 case). It is not 
immediately clear which case holds true; however, we 
have tried each possibility and find success with our ef- 
fective model if the dimension of the condensing operator 
is taken as p — A/ij. 

Now it is natural, on dimensional grounds, to include 
an additional term in the potential for our order param- 
eter </> coming from the p dimensional theory 



(60) 



AV eff (<P) = a m <f /(f- A "0 

Again aiR is assumed to be positive. 

If fp < Air < p, then the quartic contribution to 
the effective potential dominates over the term from the 
infrared theory. In this case, minimizing the potential 
yields the standard mean-field critical exponent. How- 
ever, if < A/fl < |p, then we find the non-mean-field 
exponent 



(O - O c f = (O - O c ) 



P- A IR 



(61) 



When the bound Air < A c = 3/4p is satisfied, the con- 
densate scales with a non mean-field exponent. 

For the example of the magnetic field competing with 
our operator O in the D3/D5 system. We can find the 
expected critical exponent 
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We have shown that in the D3/D5 system with a mag- 
netic field, density and a phenomenological operator O 



5 This result is based on unpublished work by KJ with T. Faulkner. 



log 



if 



f 2 




of Air and found the same basic result. The scaling of 
the free energy is reproduced by the effective potential 
Eq. (63) for the corresponding exponent j3. This shows 



the strength of the effective potential analysis which does 
not need to know even the dimension of the operator O. 



The effective potential Eq. ( 60 ) can be generalized for 



V d 



cases where the emergent theory has different scaling 
symmetries than in the deformed D3/D5 system. The 
structure is essentially the same: a Landau- Ginzburg po- 
tential analytic in 4> 2 and a non-analytic piece stemming 
from the infrared theory. In the next section we will show 
such an application to the D3/D7 system. 



FIG. 7: A log-log plot of the difference of free energy, AF, in 
the broken phase at zero temperature as a function of O. The 
dots indicate numerical data at magnetic field B — 0.95B C 
and the choice A = 1, and the line a numerical fit. The free 



energy scales as AF ~ — (O — O c 



in the broken phase. 



where we set p = 1 and the B range comes from 1/2 < 
A IR < 3/4. As a result 1/2 < j3 < oo. For B < J^d, 



13 = 1/2. We have plotted the result Eq. <|62J) over our 
numerical results in Fig[2]and they match the numerical 
results perfectly. 

In many ways Eq. ( |60[ ) is the primary result of this 
work. The non-trivial emergent theory leaves a finger- 
print on the effective potential which can lead to rich 
phase diagrams even in the strict N — > oo limit. For 
example, if Air < 5/6 then a non mean-field tri-critical 
point is realized by driving either ay (for Air > 3/4) 
or aiR (for Air < 3/4) negative while keeping the other 
positive. In the first case, the terminating line of second- 
order transitions has mean-field exponents while in the 
second it does not. 

We can further test our effective potential by measur- 
ing the free energy near the transitions we identified in 
Sec. |III A 3| if we write the effective potential as 

V eff (<i>) =V + a 2 (O c - 0)4> 2 + W 2+ ?, (63) 



In the broken phase this gives 

4> ~ (O - O c f , AF ~ (O - O c ) 1+2f3 , 



(64) 



giving a prediction for how the free energy should scale 
across the transition. 

Recall that in order to measure the free energy, we 
compute (minus) the regularized bulk action. We do this 
numerically, employing the methods of [51 . We plot 
some representative results at relatively large magnetic 
field, B = 0.95B C and the choice A = 1 in Fig. [7} Nu- 
merically, we find that the condensate scales with expo- 
nent /3 = 1.20 and the free energy in the broken phase 
as AF/Af — -(O - O c ) 3 38 . which indeed reproduces 
the scaling of the free energy for j3 — 1.20 to within 
1%. We repeated this analysis for many different values 



V. BKT AND NON MEAN FIELD 
TRANSITIONS IN THE D3/D7 SYSTEM 

As another example of our phenomenological analysis 
and effective theory methods let us finally return briefly 
to the D3/D7 system in Eq.(19). That system has a 



magnetic field, density and a phenomenological operator 
O present. 

Looking at linearized fluctuations around the symmet- 
ric L = embedding, there are three different infrared 
limits depending on the value of A. At large A > 3, the 
field L/p fluctuates unstably in the infrared. For A < 3, 
L/p fluctuates as a stable massless scalar in AdS 2 . How- 
ever, for the particular case A = 3, L/p fluctuates as a 
m 2 = -30 2 /(J 2 + O 2 ) scalar in AdS 2 . Adjusting the ra- 
tio O/d then tunes the dimension of the scalar operator 
dual to L in the infrared theory to be 



1 



Air = 



cP-llO 2 
d?+0 2 



(65) 



This last case is perhaps the most interesting. No- 
tably, it corresponds to the case where the density and 
simulated deformation have the same dimension. As with 
the D3/D5 system at nonzero magnetic field and density, 
increasing O at fixed density will trigger a chiral holo- 
graphic BKT transition at d/O c — y/Tl as the field L/p 
violates the BF bound in the effective AdS 2 region. On 
the other hand, at smaller values of O we can presumably 
drive a second-order chiral transition with the magnetic 
field. There will be a regime of non mean field transitions 
for some intermediate values of O. The choice A = 3 is 
therefore analogous to our phenomenological D3/D5 sys- 
tem with the roles of the magnetic field and O reversed. 

We can apply our effective field theory Eq. ^ to this 
case too. It predicts the critical exponents 



I d 2 + 2 

d 2 - no 2 



- i 




-d < O < J—d 
43 V 11 



(66) 



VI. DISCUSSION 

We now summarize our results. For our phenomeno- 
logical D3/D5 setup with magnetic field, density, and a 
third control parameter we find the non-trivial phase dia- 
gram in Fig. [I] The new ingredient is that tuning a third 
control parameter can lead to a line of second-order tran- 
sitions. Moreover, the critical exponents of these transi- 
tions do not appear to depend on the details of the de- 
formation. Rather they are functions of the dimension, 
A//?, of the operator dual to the embedding function in 
an emergent infrared theory. This dimension is tuned by 
the equal-dimension (in the UV theory) control parame- 
ters density and magnetic field. 

We have measured or computed four of the critical ex- 
ponents along this line. In Sec. |III A 3[ we numerically 
measured the condensate in the broken phase and found 
agreement with an analytic function of the infrared di- 
mension as in Eq. ([!]) . In terms of the magnetic field and 
density, the exponent f3 is given by Eq. (62 1. At smaller 



B/d (below y/Z/29), j3 is simply 1/2. We also computed 
the free energy in the broken phase and found that both 
its scaling and j3 follow from an effective potential Eq. 
From this we can also compute the critical exponent 7, 
which is related to the scaling of the susceptibility, 



dm 



(o-o c 



(67) 



since the hyper mass m is conjugate to the condensate ■ 
The effective potential with a mass becomes 



V m , e ff(<j>) = V e ff((f>) + m<j). 



(68) 



The exponent 7 is computed from the effective potential 
to be the mean-field value 7 = 1. 

At nonzero temperature, however, all non mean-field 
scaling is lost and the effective potential becomes an or- 
dinary quartic polynomial. The exponents f3 and 7 are 



1/2 and 1 respectively. The temperature destroys the 
non mean-field scaling, just as in the holographic BKT 
transitions. 

We also computed the low-energy behaviour of the two- 
point function of the condensate in Sec. |III C[ Fluctua- 
tions of the condensate correspond to time and spatially- 
dependent fluctuations of the bulk field L. At small 

w <C y/B, Va, the equation of motion for 9 ~ L/p re- 
sembles that of a scalar in AdS2- The equation is solv- 
able there and we match it to the physics outside the 
infrared region. This enables us to show that at T = 
the dynamical critical exponent is non mean-field, 



2 

2Am 



1 



dependent only upon Air. Moreover it returns to the 
mean field value z — 2 at finite temperature. 

In the end, our main results are dependent upon the 
details of the emergent infrared theory. As a result we 
expect them to hold in a much wider class of problems, 
including conformal theories with three control parame- 
ters as discussed. It would be nice to test this picture in a 
purely field theoretic context. In particular, it would be 
extremely interesting to realize both the emergent the- 
ory as well as the sort of phase diagram we identify in a 
conformal theory without a holographic dual. 
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